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Abstract. Slow roll reconstruction is derived from the Hamilton- Jacobi formulation 
of inflationary dynamics. It automatically includes information from sub-leading 
terms in slow roll, and facilitates the inclusion of priors based on the duration on 
inflation. We show that at low inflationary scales the Hamilton-Jacobi equations 
simplify considerably. We provide a new classification scheme for inflationary models, 
based solely on the number of parameters needed to specify the potential, and provide 
forecasts for likely bounds on the slow roll parameters from future datasets. A minimal 
running of the spectral index, induced solely by the first two slow roll parameters (e and 
rj) appears to be effectively undetectable by realistic Cosmic Microwave Background 
experiments. However, we show that the ability to detect this signal increases with 
the lever arm in comoving wavenumber, and we conjecture that high redshift 21 cm 
data may allow tests of second order consistency conditions on infiation. Finally, we 
point out that the second order corrections to the spectral index are correlated with 
the infiationary scale, and thus the amplitude of the CMB B-mode. 
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1. Introduction 

Inflation is an elegant explanation for tlie large scale appearance of our universe. 
Causally connected regions of space are swept outside the Hubble horizon during a phase 
of accelerated expansion and cross back during a later epoch of regular expansion. The 
inflationary predictions of a flat universe, along with a Gaussian, adiabatic and nearly 
scale invariant power spectrum are perfectly consistent with current data. Moreover, the 
observed anti-correlation between temperature and polarization at large scales provides 
further support for the inflationary origin of the initial density perturbations. 

While inflation predicts the overall form of our visible universe, we have very little 
understanding of the physical mechanism that generates the accelerated expansion. 
Consequently, there is considerable interest in "reverse engineering" the inflationary 
potential from astrophysical data. The primordial power spectrum can be written 
as a function of the slow roll parameters. In their simplest form, the slow roll 
parameters are expressed as derivatives of the potential, so measuring the spectrum 
to arbitrary precision would yield a Taylor expansion of the potential. Unfortunately, 
the simplicity of this scheme is undermined by the practical challenges that arise during 
its implementation. Recall that the evidence for any scale dependence in the primordial 
spectrum is still preliminary, although there is indeed growing evidence the spectrum 
is red, or has diminishing power at short scales [Ij. Reconstructing the potential 
requires reliably distinguishing between models which all predict slightly different red 
spectra, which would need an improvement in the accuracy of cosmological parameter 
determinations by a further order of magnitude. Given both cosmic variance and the 
practical challenges of foreground subtraction it is not clear that this is possible, even in 
principle. This is a particularly pressing problem if we are limited to data derived from 
the CMB [Cosmic Microwave Background] and LSS [Large Scale Structure]. These two 
sources of data are sensitive to a range of scales that differ by a total factor of perhaps 
10^ but in simple models of inflation the primordial spectrum changes very slowly, 
undermining the power of any reconstruction program. Despite these challenges, there 
are powerful motivations for theoretical studies of reconstruction. The first is that the 
physical basis of inflation (if inflation is, in fact, the source of the primordial fluctuations) 
is one of the most important open questions in all of cosmology, and answering it is likely 
to shed light upon particle physics at very high energies. Secondly, a number of proposed 
experiments aim to measure the fundamental spectrum to very high precision; we wish 
to determine their ability to constrain the overall inflationary parameter space, even 
though many of them will take decades rather than years to implement. 

The reconstruction of the inflationary potential was was first discussed in the 
early 1990s [21 El HI [5]. This paper builds on the slow roll reconstruction algorithm, 
proposed and implemented by Easther and Peiris P, [3 E], which grew out of a Monte 
Carlo approach based on the inflationary flow equations [9l [TOl [TTl [121 IS]- Originally 
reconstruction was based on taking the measured values of the spectral indices (and their 
running), solving for the slow roll parameters, and deducing the form of the potential 
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[2]. Since we are interested in the inflationary parameter space, there is no need to 
compute the spectral indices, as these have no fundamental significance. Rather, one 
can include the slow roll variables directly in the cosmological parameter set and bound 
them using Monte Carlo Markov Chain fits to cosmological data p!5| [6l [7] . Slow roll 
reconstruction is thus an optimal approach to recovering the infiationary potential from 
data, since it makes use of all available information. Moreover, while it makes use of the 
slow roll expansion, we need never employ the slow roll approximation, as the truncated 
Hamilton- Jacobi hierarchy can be solved exactly [I6]. Slow roll reconstruction thus 
captures all the correlations between the slow roll parameters, including those at second 
and higher order 1^1 

A number of other approaches to reconstruction have been proposed. In particular. 
Leach and collaborators write the spectral indices (and their running) in terms of the 
slow roll parameters at a fixed pivot [IHl [20], El] . When used with data which probes a 
small range of scales, this approach is functionally identical to slow roll reconstruction, 
since the slow roll parameters are effectively constant. However, we will show that as 
the lever arm in wavelength becomes large, the scale dependence of these parameters 
can be significant. In principle, one could account for this running by computing the 
perturbation spectrum using expressions that include higher order corrections in slow 
roll, but this approach will become algebraically cumbersome at some point. Conversely, 
[22] explores constraints on the slow roll parameters imposed by demanding a sufficient 
period of infiation. Finally Cline and Hoi look at reconstructing infiationary models 
with significant running within the Hamilton- Jacobi formalism [23]. Since slow roll 
reconstruction implicitly predicts the form of the potential it naturally makes use of 
this information, and one can restrict fits to parameter values that allow a sufficiently 
long period of infiation. 

The aim of this paper is to forecast the parameter constraints we can expect 
from slow roll reconstruction when it is applied to future datasets, and to explore the 
differences between slow roll reconstruction and fits to the usual spectral parameters. In 
the process, we show that low scale infiation (i.e. energies significantly below 10^^ GeV) 
is described by a simplified set of Hamilton- Jacobi equations, as the first parameter 
(e) is effectively absent from the dynamical system. We present a new classification 
scheme for infiationary models based on the number of free parameters needed to 
specify the potential - rather than its shape - and show this is naturally related to 
the (truncated) slow roll hierarchy. We then use Fisher matrix calculations [211 [25] 
to explore how constraints on the infiationary parameter space improve with the lever 
arm in wavelength probed by the dataset|§| With a purely Gaussian likelihood function 

I As explicitly implemented in [6l [7], slow roll reconstruction uses approximate expressions for the 
perturbation spectrum, but this is simply a matter of convenience 0. One can always solve the 
perturbation mode equations numerically, and this approach has been explored [171ll8j. As we will see, 
computing the spectrum with the scale dependent slow roll parameters yields a very good match to the 
exact calculation. Numerically evaluated spectra eliminate a source of uncertainty in the parameter 
constraints, but are not justified by the quality of presently available data. 

§ We utilize CAMB, http://camb.info/, [26^ for calculations of the CMB spectra and the matter power 
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(which is an assumption of the Fisher matrix formahsm) the CMB alone does not 
probe a large enough range of wavelengths to produce results that differ significantly 
from those obtained via fits to the standard spectral parameters. However, if one has 
access to the primordial spectrum at very short scales (e.g. via high redshift 21 cm 
measurements), terms that are second order in slow roll may become significant. In 
practice, the likelihood is far from Gaussian, so this analysis is effectively a worst case 
scenario. In particular, slow roll reconstruction allows us to include information about 
the duration of inflation into parameter estimates; these constraints are not captured by 
a Fisher matrix analysis. Models with a significant tensor spectrum generically require 
a longer period of inflation than those with an unobservable tensor component; we show 
how this information - along with a very mild prior on the post-inflationary equation 
of state - can give further leverage to slow roll reconstruction. 

Since we are considering constraints on the inflationary parameter space, we 
are primarily interested in models for which f^Tot = 1; the flatness of the spatial 
hypersurfaces is a key prediction of inflation. This is an implicit assumption of slow 
roll reconstruction, since the Hamilton- Jacobi equations are derived after ignoring the 
spatial curvature term in the Einstein equations. If i^Tot differs slightly from unity 
today the longest modes will have left the horizon just as inflation began. Consequently, 
slow roll reconstruction (as currently implemented) is only self-consistent when inflation 
lasts long enough to suppress any transients associated with the pre-inflationary initial 
conditions before modes which contribute to the quadrupole leave the horizon. We also 
assume that the entire primordial perturbation spectrum was generated during inflation, 
and contains no significant contribution from cosmic strings, or other non-inflationary 
mechanisms. Finally, we treat the dark energy as a pure cosmological constant, although 
relaxing this assumption would not significantly modify our key conclusions. 

The CMB temperature peak morphology is very well understood [2811291, [30 t [311 [32] . 
along with the importance of the E and B polarization modes. The ith multipole of 
the CMB anisotropy corresponds to a wavenumber ^ 2/{Hq€), so the CMB probes 
scales from ~ 2 x lQ~^h Mpc"^ at the quadrupole to A; ~ O.lS/i Mpc"^ at ^ ~ 1500. 
Beyond £ ~ 1500, the primordial Ci decay sharply due to Silk damping, which reflects 
the finite width of the surface of last scattering. It is thus difficult to measure the 
primordial CMB at shorter scales. Moreover, foregrounds and secondary anisotropics 
typically grow in amplitude at smaller scales. In what follows, we will consider perfect 
CMB measurements out to ^max of up to 2,500 in order to explore the lever arm yielded 
by data over a large range of angular scales. However, given the difficulties associated 
with foreground subtraction, these calculations are essentially gedanken experiments. 

On the other hand, we are not limited to CMB information alone. The information 
contained in LSS data is largely orthogonal to that in the CMB, breaking many 
parameter degeneracies. This 'cosmic complementarity' is well known [331 [Ml [271 [35] , 
and the slow roll reconstruction can obviously make use of this data. While the 

spectrum, and the two sided derivative methodology outhned in [27j . 
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linear regime of structure formation is very well understood, recovering the primordial 
spectrum at length scales which have undergone nonlinear evolution is a challenging 
task. The smallest scale still in the linear regime probed by low low-redshift LSS 
experiments is not wildly different from that probed by high resolution measurements of 
the CMB. In what follows, we are primarily interested in the differences between slow 
roll reconstruction and analyses based on the usual spectral variables, and adding LSS 
data would improve both approaches. Furthermore, forecasts of parameter uncertainties 
for LSS data depend significantly on the design of the experiment. Consequently, in 
this analysis we have focussed on CMB data alone. One can probe smaller scales by 
looking at very high redshift data, when these modes were still in the linear regime. 
Options in this area include very deep galaxy surveys, or Lyman-a experiments. More 
speculatively - but with far greater potential power - the high redshift 21 cm background 
may allow probes of the primordial spectrum at very small wavelengths [36]. We plan 
to look at constraints derived from combined fits to CMB and 21 cm data in a future 
paper, but the calculations here suggest that such a measurement might permit tests 
of higher order inflationary consistency conditions [371 EH HI]- Finally, a BBO style 
experiment is sensitive to the primordial tensor spectrum at solar system scales, and 
slow roll reconstruction could take advantage of this data to put exquisitely accurate 
constraints on the inflationary potential. 

This paper is organized as follows. In ^ we review inflationary dynamics, the slow 
roll approximation and our construction of the primordial spectrum. In ^ we consider 
analytic models with low inflationary scales, and construct a new classification scheme 
for inflationary models in ^ In ^ we review the Fisher matrix formalism for CMB 
anisotropy experiments, and describe the resulting forecasts in ^ We conclude in ^ 
Finally, in an Appendix we give a more detailed analysis of the slow roll dynamics at 
low inflationary scales. 



2. Inflation and the Primordial Power Spectrum 
2.1. The Background 

We use the Hamilton- Jacobi formulation of inflationary dynamics, expressing the Hubble 
parameter as a function of 0, rather than as a function of time. Thus H = H{(j)), and 
we assume that is monotonic. The equations of motion are [HI [39l HHl HH HJl 133] 



4:71 



H\<P). (1) 



[H\<P)f-^^H\ct.) = -^-^V{ct.). (2) 



Primes denote derivatives with respect to the field, while an overdot denotes derivatives 
with respect to coordinate time. Equation ([2]) is the Hamilton- Jacobi equation, and 
describes inflation in terms of the Hubble parameter, H{(j)), rather than the potential, 
V{(j)). The Hubble parameter, being a geometric quantity, describes the spacetime 
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dynamics, whereas particle physics constructions predict V{(f)). We can thus discuss 
slow roll inflation without specifying the particle physics that generates inflation. The 
HSR [Hubble slow roll] parameters ^Xh are defined by the infinite hierarchy of differential 
equations [10]: 
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> 1. 



(3) 
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The usual slow roll parameters are t] = ^Xh and ^ = "^Xh- If we truncate the hierarchy, 
so that ^Xh = for all £ > M at some 0O; then these ^Xh vanish everywhere. When 
truncated at order M, the hierarchy can be solved exactly to obtain [16] 
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where the -B„ are specified by the initial values of the HSR parameters: 
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The subscript refers to their value at the moment the fiducial mode fco leaves the 
horizon (when aH = ko) and = 0o = 0. In this analysis we set ko = 0.05 Mpc~^, 
which corresponds to £ ~ 500 in the CMB. 

Substituting equation ([3]) into equation ([2]) gives the potential 
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while N , the number of e-folds before the end of infiation is 
dN in H 1 
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The scale dependence of the slow roll parameters follows from equations ([3]) and ([8 
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(9) 

|: 

(10) 

(11) 
(12) 



and the truncation property mentioned above can be derived from the last equation. 
Note that these formulae use N as their independent variable. 

Once we have specified the values of the slow roll parameters at some fiducial scale, 
we may "fiow" to any other scale by using the slow roll hierarchy. The scale dependence 
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of the slow roll parameters ensures that any non-trivial correlation between the slow 
roll parameters will only strictly apply at a single value of (p. For instance, if might 
assume the existence of a special point where e and t] were non-zero, while 'A^ = for 
i = 3, ■ ■ ■ , N — 1 and ^ Xh 7^ (see ^1] for an example). In principle, one could analyze 
this type of model by adding ko to the parameter set in order to marginalize over the 
model-dependence in the mapping between k and (p. 



2.2. The Perturbations 



The scalar and tensor perturbations obey 



z dr'^ J 

I d^z 
z dr"^ 



(13) 



(14) 



Here r is conformal time, Uk are the Fourier modes of the gauge invariant Mukhanov 
potential describing the intrinsic curvature perturbation, Vk are the analogous Fourier 
modes for the tensor perturbations, while z = a(f)/H for scalar perturbations, and z = a 
for tensor perturbations. The power spectra are 

(15) 
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A first order expansion about tlie exact solution for power law inflation gives [IHl El [7] 



[1 - (2C + l)e + Cr/]2 / H 



Vh = [l-(C+l)e]' 
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where C = —2 + ln2 + 7 ^ —0.729 and 7 is the Euler-Mascheroni constant. The power 
spectrum is normalized at k^ by setting 

[1 - (2C + l)eo + Criof f 



(19) 



where eo, ?7o and Hq are the values of these parameters at /cq- 

One can always solve the perturbation evolution equations numerically |17l H5| HSJ 
50] . This approach is implemented in p!71 [18] who parametrize H^cf)) as a polynomial 
of finite order - which is identical to slow roll reconstructionJUj since the truncated flow 



II Note that [TTl [18] also consider fits to a polynomial expression for Hamann, Lesgourgues, 

and Valkenburg [18] present constraints on the first three slow roll parameters from WMAP3 [5T] and 
ACBAR [52] derived with two different approximations to the spectrum, as well as the numerically 
evaluated mode equations. Their analysis shows that the difference between these results is entirely 
explained by the implicit priors on the duration of inflation, and not the accuracy with which P{k) is 
evaluated. The ability to include constraints based on the duration of inflation is a key feature of slow 
roll reconstruction, and we return to this topic below. 
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Figure 1. The top panel shows the primordial spectra for two parameter choices; 
e = 0.01, 77 = 0.01, C = (upper) and e = 0.01, tj = 0.01, ^ = 0.001 (lower, curved). 
The bottom panel shows the difference between the spectra, the dashed curve applies 
to ^ = 0.0 while the solid curve corresponds to ^ = 0.001. For extreme values of ^ these 
discrepancies can become large, but these models typically have a very low number of 
e-folds. 

hierarchy is solved exactly by a polynomial in H{(f)) \T5]. Since equation (fTTll uses the 
scale dependent slow roll parameters - which are matched to a value of k by solving 
equation ([9]) - it accurately tracks the exact spectrum for the combinations of slow 
roll parameters one is likely to encounter in practice. We plot the difference between 
equation (fT7|) and the exact spectrum for two representative sets of slow roll parameters 
in Figure [H In this paper, we perform our calculations using equations (fTTI) and (fT8|) . 
but a numerically computed spectrum would banish the last vestiges of the slow roll 
approximation from this analysis, and might be justified when dealing with very high 
quality data, especially if the running in the spectral index turns out to be non-trivial. 

From (fT7|) and (fT8|) . one can also recover expressions for the usual spectral indices 
and their scale dependence jl4j : 

Us = 1 + 2r/ - 4e - 2(1 + C)e' - ^(3 - 5C)er] + ^(3 - C)^, (20) 
r = 16e[l + 2C{e-ri)], (21) 
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where C = 4(7 + In 2) — 5 and is not to be confused with C above. The expressions are 
not used during slow roll reconstruction, but are employed when we make comparisons 
to the empirical characterization. 



2.3. The Duration of Inflation 

We need make no explicit assumption about the duration of inflation. Other HSR 
based analyses often pick a specific set of parameters, then integrate forward in time to 
find the moment at which inflation ends, and then work backwards from this point to 
find the slow roll parameters some fixed number of e-folds before the end of inflation 
frU\ |TT[ |12l [13] . By running Markov Chains with the slow roll variables in the parameter 
set we learn their values at the moment when cosmological scales were leaving the 
horizon. We can then evolve them forwards in time, to discover the point (if any) at 
which inflation ends, which occurs when e = 1. In three parameter fits to current data, 
one finds considerable support for a running index, which requires large values of ^, 
leading to models with ~ 15 or even less |71 E]. If we assume that we have enough 
slow roll parameters to describe the inflationary potential this situation is clearly not 
self-consistent. One can posit the existence of a secondary period of inflation, or expand 
the parameter set in order to make ^ strongly scale-dependent, but either of these 
"solutions" implies that the set {e, r/,^} is incomplete. Conversely, if we find a large or 
even unbounded number of e-folds for a given set of slow roll parameters, the implication 
is that inflation ends via some sort of hybrid or waterfall transition [531 EU > where the 
field evolves in a direction orthogonal to the original inflaton trajectory. These two 
scenarios are analogous to following a path that leads down a mountainside to a valley 
or walking along a gentle plain and then encountering a sharp cliff, respectively. In 
the former case, the end of inflation can be extrapolated from the slope of the potential 
some distance away from the minimum, while in the latter case the inflaton field receives 
little or no advance warning that inflation is about to end. 

The primordial tensor amplitude is set by the energy scale of inflation. Assuming 
slow roll, the scalar amplitude is given by [55] 

Vs, = l^i. (24) 

where Vs,o — 2.6 x 10^^ [20j and Vk is the value of the potential as the mode k leaves 
the horizon. This energy scale can be translated into a constraint on the number of 
e-folds needed to reproduce the observed universe [15] 
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where Vend is the energy scale at the end of inflation, p^e/i is the energy density at 
reheating and ao-ffo is the expansion rate today. As pointed out by Kinney and 
Riotto [56j, the undetermined parameters in A^(A;) induce a theoretical uncertainty in 
associating inflationary parameters with a potential. 

By writing down equation (l25l) we have assumed that the universe is matter 
dominated after the end of inflation and then radiation dominated between the end 
of reheating and matter-radiation equality, which amount to assuming that the effective 
equation of state parameter varies between w = and w = 1/3 between the of inflation 
and nucleosynthesis. If we include matter with w > 1/3, such as a kination field |57j . 
N{k) will increase. The minimum possible value of in a decelerating universe is —1/3, 
and if —1/3 < w < 0, N{k) can be substantially reduced. Note that the term involving 
Vk implicitly contains e. However, in slow roll inflation Vk is not expected to differ 
considerably from Vend, and thus this term is not considered to be important. 

Slow roll reconstruction allows one to include priors based on the duration of 
inflation when fitting to the slow roll parameters. In particular, the duration of inflation 
depends very strongly on the value of C, [8J, and even the modest requirement that 
> 30 induces stringent bounds on the parameter space [6]. However, studying 
equations and (125!) shows that the value of depends strongly on e - while we 
can certainly have inflationary models with A^ ~ 30, this requires a tiny inflationary 
scale, and thus a miniscule value of e. Conversely, if e ~ 0.01, the inflationary scale 
must be comparatively large, and A^ ~ 50. Consequently, we can implement slow 
roll reconstruction with an e dependent constraint on the number of e-folds by using 
equation to write as a function of e and the amplitude of the primordial 
perturbation spectrum. This requires a mild prior on the post-inflationary equation 
of state - for instance that w > after the end of inflation. This implicitly rules out 
exotic scenarios such as a phase in which the universe is dominated by a frustrated 
network of cosmic strings [58], but does not strongly constrain the post-inflationary 
universe. Finally, a principal lesson of [18j is that all implementations of slow roll 
reconstruction make implicit assumptions about the minimal permissible duration of 
inflation. Consequently, slow roll reconstruction is most transparent when an explicit 
prior on A^ is included in the chains, even if that bound is very mild, and we will pursue 
this topic in a separate paper. 

3. Slow Roll in the Low-e limit 

A primordial gravitational wave spectrum is often described as the "smoking gun" of 
inflation, since there is no credible alternative mechanism for generating long wavelength 
gravitational waves. Detecting this signal would determine the overall energy scale of 
inflation since Vh^ H'^. Intri guingly, string theoretic inflationary scenarios generically 
predict a negligible amount of primordial gravitational radiation [59], or r < 10~^°. In 
this case, e must also be tiny, since r ~ 16e. Conversely, models with algebraically 
simple potentials typically have r > 0.001 [60]. The constraint on stringy models can 
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be understood in terms of the Lyth bound [SI], which requires that the inflaton's total 
excursion be sub-Planckian to ensure that its potential is not dominated by contributions 
from higher-order operators. 

The current observational bound on r is relatively weak, r < 0.3 ^51j, but upcoming 
CMB experiments may push this value down to r ~ 0.01 [62l ESl Elj. However, if r is 
many orders of magnitude smaller than unity it may be unobservable by any conceivable 
experiment [65]. While e must be tiny in these models, there is no equivalent constraint 
on 1], and in many supergravity or stringy models the challenge is to ensure that \ri\ < 1. 
Consequently, we assume 1 |?7| ^ e, and write the flow equations as 

^ = (e~i)n'x„ + '»XH. (28) 

Here ^ drives the evolution of t] and ^Xh drives the evolution of ^ (if we could measure 
it), and one can still truncate the hierarchy. We can also estimate e, 

^ = 2er^ + 0{e% (29) 



which can be solved to give 

e(iV) = e(iVo)exp 



N 



(30) 



7]{N)dN . 

'No 

Since -C 1, if e is initially small, it stays small and then rises super-exponentially 
as 7] ~ 1. Conversely, if 1] is also vanishingly small, we cannot usefully extend this 
truncation, since e is unique in that both of its source terms in the slow roll hierarchy 
are suppressed when it is small. To illustrate, suppose that the first m — 1 slow roll 
parameters are close to zero, but "^Xh 7^ 0. However, "^Xh generates "Aj^ 7^ for all 
n < m, so this condition amounts to choosing the initial conditions at a very special 
point in the potential. The coupling between "^Xh and "^~^Xh means that after one 
e-fold "^~^Xh will be at least as large as "^Xh- One may proceed inductively to see that 
after m e-foldings (and probably sooner) all of the slow roll parameters except for e will 
differ significantly from zero. 

When we drop e, we can solve the slow roll hierarchy exactly if it is truncated at a 
relatively low order. These solutions are explored in Appendix A Consider the specific 
case ^ = 0, ?7 <^ 1 and e -C 77, so dr]/dN ^ and ri{N) ^ 77(0) = 770, which corresponds 
to the "Low-e 1-Parameter" model defined in the next section. From equation (l22l) it 
follows that Us is effectively constant over cosmological scales. Equation (JTOl) has the 
solution 

eo + (^70 - eo)e ^''"^^ 
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Figure 2. Slow roll hierarchy truncated at ^ = 0. The plot shows the evolution of e, 
the red curve is the exact result, the dashed blue curve is the approximation with rj 
constant. We have set r] = —0.02 and e = 10~^. 



where we have set = at the fiducial scale for convenience. In this approximation, 
we can solve exactly for where infiation ends: 



and we recall that unless rjo < 0, infiation will continue indefinitely in this scenario. 

If we have e ~ 10~^°, and r] = —0.02 (e.g. Ug ^ 0.96), we have infiation ending 
~ 500 e-folds after the fiducial scale has left the horizon while e ~ 10~^ gives a futher 
~ 190. As illustrated by Figure [2], the approximate solution is typically good up 
until the last couple of e-folds and leads to a slight overestimate of A^. Consequently, 
if 1] ^ —0.02, and e is very small, infiation either ends via a hybrid transition when 
it encounters an abrupt cliff in the potential, or at least one higher order slow roll 
parameter has a non-trivial value. Further, looking at equation (122!) . if e is tiny and 
^Xh = for £ > 2 the running is entirely dominated by ^ and is thus negligible if ^ = 0. 
Conversely, if e ~ 0.01, some running will be generated by the term in the slow roll 
hierarchy, even if ^ = 0. 

We also see that ^ quickly dominates the dynamics when it is of the same order 
of magnitude as t], as happens when we try to match the central value of the running 
derived from WMAP data. As noted by [8] the effect of this is to reduce the number of 
e-folds to an unacceptably small level. If the running is significant, then we would need 
at least one further slow- roll parameter in order to ensure that ^ itself is scale dependent. 
Unless such a scenario is a prediction of a well-motivated model, this situation would 
justifiably be regarded fine-tuning. 




(32) 
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Class 


HSR 


Spectral 


Shorthand 


High-e 1-Parameter 


e 


ris — 1 = tT = —lut, a ~ 




Low-e 1-Parameter 


V 


7^ 1, r 0, a 


ACDM 


High-e 2-Parameter 


e, 7] 


- 1 7^ |r = -|nt, « ^ 


ACDM+r 


Low e 2-Parameter 


V, ^ 


ris 7^ 1, r «i 0, a 7^ 


ACDM+a 


High-e 3-Parameter 




ns - 1 = |r = -\nt, a 7^ 


ACDM+r+a 


Low-e 3-Parameter 




7^ 1, r ^ 0, a 7^ ... 





Table 1. We show the non-trivial variables for models with up to three HSR 
parameters, and the equivalent set of variables written in terms of the scalar and tensor 
indices. The final column denotes the corresponding shorthand for the resulting model. 
The canonical ACDM case corresponds to assuming that e ^ |?7|. In this column r 
refers to adding a tensor spectrum, while a denotes the running of the spectral index. 

4. Parameter Counting and Model Classification 

We now describe a classification scheme for infiationary models, based on the slow 
roll hierarchy. In the past, it has been tacitly assumed that e is the one parameter 
that cannot be dropped from the slow roll hierarchy. Consequently, all previous 
implementations of slow roll reconstruction have allowed for a tensor component in the 
CMB, whereas the minimal ACDM parameter set includes only scalar perturbations. In 
the language of slow roll, the ACDM parameter set is thus equivalent to assuming that 
e < 10~^, in which case e plays no direct role in determining the observable properties 
of the present universe, and we need to include it in our chains. We thus propose the 
following scheme: 

• Low-e, A^-parameter. The tensor-scalar ratio is assumed to be immeasurably small, 
and we truncate the slow roll hierarchy at ^+^A. 

• High-e, A^-parameter. The tensor-scalar ratio can be measurably different from 
zero, and we truncate the slow roll hierarchy at ^A. 

This is distinct from the "large field / small field / hybrid" zoology (see e.g. |12j). 
Models with a large variation in necessarily have a non-trivial e [61] [66] , whereas 
the small field and hybrid cases are distinguished by the sign of rj. However, when ^ is 
nontrivial, r] can change sign during the course of infiation, so "shape-based" taxonomies 
are best avoided with general potentials. Note that this classification does not specify 
how infiation ends, whether via a hybrid transition, or a violation of slow roll - the 
valley or the cliff. 

Within this schema, we can consider A^ parameter models. A Low-e 1-Parameter 
model is almost entirely equivalent to standard ACDM, since 7] is the only nontrivial 
slow roll parameter. In this case the running in r] is effectively zero, which leads to a 
constant spectral index, with no tensor component. Conversely, a High-e 1-Parameter 
model has no obvious analog, since ra^, r, Ut and a are all specified in terms of a single 
parameter. Physically, setting = means that -f^(0) is a pure quadratic function. 
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Since e increases with time, a very small value of e at CMB scales would imply that 
inflation must have ended via a sudden hybrid-style transition. 

With two parameters, a High-e model corresponds to ACDM+r, where the tensor 
spectrum obviously has an inflationary prior that links the amplitude and spectral index. 
Moreover, as we will see later, the flow equations imply a weak scale dependence in e 
and rj. Thus, the running of the effective spectral index is non-zero - and for very high 
quality data, this increases the leverage we can obtain from slow roll reconstruction. 
Conversely, a Low-e, two parameter model corresponds to ACDM+a, but we now must 
be careful to ensure that for comparatively large and positive values of ^ the total 
duration inflation is self-consistent. If we do wish to consider a large running at CMB 
scales, our one recourse is to add ^Xh to our parameter set, leading to either a High-e 
4-Parameter model, or a Low-e 3-Parameter model. In this case the tilt is transient, but 
we are left with a potential that has several free parameters. The classification scheme 
is summarized in Table |H 

As we noted above, this analysis does not specify the mechanism that ends inflation 
- that is, whether we are rolling toward a cliff or a valley. However, since slow roll 
reconstruction effectively specifies the overall form of the potential, we can add further 
cuts to our parameter space by insisting on a sufficient overall duration of inflation, 
as described at the end of Section 2. This is not possible when using the spectral 
indices and amplitudes, since these variables contain no information about the duration 
of inflation. 



5. Error Estimation and Fisher Information 



For our error forecasts we make use of the Fisher information matrix [23] , a measure of 
the width and shape of the likelihood function around its maximum. 



where L = In £ and the denote model parameters. The Cramer- Rao inequality then 
says that the minimum possible standard deviation on a single parameter, a^, estimated 
from the data is l/y/TYi- This minimum standard deviation rises to l/^{F-^)u, if all 
parameters are estimated from the same data. Previous treatments of this topic include 
[671 [62] , and we follow the formalism laid out in [62] . 

We restrict our attention to CMB data - in this paper, our principal concern is 
to compare slow roll reconstruction to fits to the spectral variables. Adding more 
information further constrains the free parameters, and thus accentuates the advantages 
the HSR formulation. Moreover, we are not considering the impact of priors based on 
the duration of inflation (which puts sharp cuts on the allowed region of parameter 
space), so what follows is essentially a worst case analysis for slow roll reconstruction. 

Observations of the CMB measure the polarization and the anisotropy of the 
temperature of the radiation in terms of spherical harmonics, from which we obtain the 
Ce for each of the spectra Cti, Cei, and Cse and the cross- correlation Cce- Assuming 
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Experiment 



fsky Frequency 6beam ctt o"p 
(GHz) (') {f,K) {f,K) 



Planck (Plk) 



143 8.0 5.2 10.8 

217 5.5 11.7 24.3 



Ideal Sat (Ideal). 



0.8 



30 - 200 



8.0 



2.2 



2.2 



Table 2. Experimental specifications for CMB satellites. 



that the CMB multipoles are Gaussian distributed and letting a denote our vector of 
parameters, with a the fiducial values, the Fisher matrix for a temperature/polarization 
measurement can be written 



We assume a gaussian beam profile and assume that foregrounds are perfectly 
subtracted. The elements of the symmetric matrix Cf^ are enumerated in |62j . 
We specify a proposed experiment in terms of at = ©fwhm/VS In 2, the Gaussian 
beamwidth, where Gfwhm denotes the "full width at half maximum" power of the 
beam. The noise per multipole is uq = a^^^Qpi^, where Qpi^ = ©fwhm = 4:71 f sky/ Npi^ is 
the beam solid angle, Npix is the number of pixels (independent beams) in the survey 
region and fsky is the fraction of the sky covered by observations [62|. The variance 
per pixel is a^^^ which can be obtained from the detector sensitivity as api^ = s/\/Nt, 
where N is the number of detectors and t the integration time per pixel. 
In this work we consider three cases: 

(i) A cosmic variance limited survey - uq = 0, (C.V.), 

(ii) The projected errors from the Planck satellite, (Plk), 

(iii) The ideal satellite experiment of ^62j (Ideal). 

For all the above experiments we take imax = 1500, unless stated otherwise. For the 
Planck satellite, we follow [27] and [67j and use only two of the ten available channels 
(143 and 217 GHz), assuming that the remaining channels have been used for foreground 
subtraction. This simple estimate leads to somewhat optimistic predictions for Planck's 
sensitivity to a primordial B-mode. However, since we are interested in the differences 
between the standard spectral analysis and slow roll reconstruction, it suffices for our 
needs. For the full sky "Ideal" experiments, we assume five identical channels at 
frequencies 30, 50, 70, 100 and 200 GHz. The satellite noise specifications are shown 
in Table [21 In all cases (except where explicitly specified) we include the effects of 
gravitational lensing of the E-mode polarization into 5-mode polarization [3T], which 
is computed within CAMB. We do not consider cases in which the B-mode has been 
delensed. We also tested our code by checking we could recover the results of [62] in the 
cases which overlap with our assumptions here. 




(34) 
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Parameterizations 

The Fisher matrix formahsm forecasts the hkely error elhpse for any given fiducial model 
- and the size of this ellipse is a function of the chosen central parameter values. For 
concreteness, we assume that the non-inflationary parameters are well described by a 
ACDM cosmology. In what follows we set the central values of the baryon fraction, Qb, 
the cold dark matter fraction, flcdm] the reduced Hubble parameter, h; and the optical 
depth to reionization, r; to their central values found in the WMAP3 concordance 
cosmology. As noted earlier, consistency with an inflationary prior requires f^Totai = 1- 
Specifically we fix f^A = 1 — ^cdm — and assume the dark energy has the equation of 
state w = —1. We ignore the neutrino mass. 

We now compare models parameterized in terms of spectral variables, (the 
amplitude, A{ko); the scalar index, n^; and the running, a; as well as the corresponding 
tensor variables) to an HSR analysis. We define our amplitude variable A(fco) via 

P{ko) = 2.95 X 10-^A{ko). (35) 

In the analysis below we always impose the simple inflationary prior on the tensor 
spectrum Ut = —r/8, where r is the ratio of the tensor and scalar amplitudes. This 
gives a considerable advantage to the spectral variables, since a generic treatment of 
the tensor modes would not assume this correlation, whereas it is implicit in the HSR 
formalism. However, it is simple to apply and used in most analyses of CMB data; so 
we adopt it here. 

The likely parameter constraints on spectral variables are well understood (see [62] 
for instance), but this is the first forecast of the expected errors for HSR parameters. We 
will show that with very high quality data the HSR formalism results in a more sharply 
peaked likelihood surface, as it can make use of higher order correlations between slow 
roll variables. Since we can place cuts on the HSR parameter space using an "e-folds" 
prior, which is not included in our Fisher matrix analysis, the following treatment 
will necessarily underestimate the strength of the slow roll reconstruction formalism. 
Ironically, this is particularly noticeable with less precise data - if the likelihood contours 
are relatively large, they are more likely to extend into regions excluded by an e-folds 
prior. 

The Lever- arm effect 

It is intuitively clear that extending the range of wavenumbers (typically measured in 
terms of the number of decades of k spanned by the data) over which we have CMB 
data will tighten the parameter bounds. However, if the spectral indices are very well 
approximated by the lowest order slow roll expressions [ug = 1 + 2r] — 4e and r = 16e), 
and e and rj do not change significantly as observable modes leave the horizon, the slow 
roll formalism amounts to a linear transformation of the spectral variables. In terms of 
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the spectral indices and running, the spectra take the form 

^ \ ns{k)-l 



where 



Pn{k) = ^ ) ' (36) 



ns{k) = Usiko) + ^ In ( ) + . . . , (37) 



and 



2 V^o 

PH{k)^Alr{£j\ (38) 

where Ut is assumed to be scale invariant. 

The extra information exploited by slow roll reconstruction comes from the 
hierarchy of consistency relations that exist between the spectral parameters at all 
orders in slow roll. Consider, for example, the case in which we have two spectral 
parameters, {ns,r}, which is analogous to the High-e, 2-Parameter inflationary model, 
specified via {e, 77}. In the spectral parametrization, a = 0, but slow roll requires 

a « -^-^ [8e^ - lOer/] -e [Be - lOr/] (1 + e). (39) 

The corresponding spectral model {ris^r} has only the term linear in Xnik/ko) in 
ln(P7^(A;)), whereas the slow roll parametrization has additional terms of the form 
(9(e, ?7, . . .) ln(A;/A';o)^. Even though the co-efficient is small, when the lever arm (in 
In/c/fco) is large the quadratic terms will make a substantial contribution to the 
spectrum. 

One can show that, to a good approximation, derivatives of the Ci with respect to 
the parameters that specify the power spectrum can be written as 
1 _ 1 dP{k) 

~/ — ^ 

Ci dai P{k) dai 

where ctj = nsidrin/dhik^e^r] etc and P{k) is the appropriate spectrum (scalar for 
temperature and E'-mode, tensor for S-mode) . Here the quantity 770 — 77* denotes the 
conformal time interval between today (770) and last scattering (77^,) and is not to be 
confused with the slow roll parameter 7/. The elements of F^j for which ctj and aj are 
both spectral parameters are approximately 

(2£ + l) d\nPn{k)d\nPT^{k) 



(40) 



^max 



2 dai daj 



(41) 



If we are working with {us, a, r} we have 



dlnCj ^dlnP^ik) ^^^r^ (42) 



dug dug \ko 

dlnCf ^dlnPnik) ^^3^ 



da da 2 \ko 

nC 
dr 



dlnCj 



(44) 
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where the last hne signifies we are neglecting the contribution of tensor modes to the 
temperature anisotropies. Assuming the consistency relation, 

0, (45) 
0, (46) 



9 In Cf 


din Ph{k) 


dus 


dus 


dlnCf 


d In Pf^ik) 


da 


da 


dlnCf 


d\nPh{k) 


dr 


dr 



The corresponding slow roll parameter set is {e, Keeping the leading terms in slow 

roll multiplying each power of In/c/fco we find 

dlnCi dlnPnik) , , \, , , / A; 

« = (-4 - 4(1 + C). - -(3 - 5C)„) In (- 



1 



lOi, - 16e + i(5 - 7C)? 



dr\ dr] 2 \/co 

1 „ 1 A , / A; \^ 2 ^/ k 



(9^ 2' ' \koJ 2! V 2 



1 , \ , f kV 1 , / ^ ^ 



where in all instances, k is understood to be k = i / {rjo — rj^,) . Similar expressions exist 
for the tensor spectrum. As noted above, for a given number of parameters, the slow 
roll parametrization adds an extra power of ln{k/ko), which is an expression of a higher 
order consistency condition [371 EHl [H] • This extra power provides more leverage on the 
likelihood space as one adds information further from the fiducial scale. 

To compare the constraining power of the two parameterizations, we consider the 
volume of parameter space contained within the Fisher ellipses, as a function of imax, 
in analogy to the metric developed by the Dark Energy Task Force [68]. For three 
slow roll or spectral parameters, this is the volume of an ellipsoid. In order to make a 
comparison between the two different parameterizations, we normalize the volume by 
its value when i^ax is set to the pivot scale which, for ko = 0.05 Mpc~^, is £ ~ 500. 
The constraining power of the parametrization determined by how quickly the volume 
contracts as information is added, and we plot a Figure of Merit, defined as the inverse 
of the normalized n-dimensional volume enclosed within the error ellipsoids. In reality, 
the CMB probes a relatively small range of scales (the largest value of In (k/ko) is likely 
on the order of 6), which limits the maximal extent of our lever arm. Moreover, the 
highest order term in ln{k/ko) in the derivatives of the Ci is multiplied by a slow roll 
parameter, which is necessarily considerably smaller than unity. Consequently, if we 
restrict attention solely to CMB data, we will need a very accurate measurement of 
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500 1000 1500 2000 2500 

MultiPole 

Figure 3. Model with e = 0.015, tj = -0.02, ^ = 0.0. The figures of merit for 
the {e,?7,$} (red) and {r, ns,a} parameterizations (black), normalized ai £ — 500, 
calculated using the full Fisher matrix with cosmic variance only. As the lever arm in i 
becomes longer, the extra information in the slow roll formalism leads to a more sharply 
peaked likelihood function and thus tighter constraints on the model parameters. 

the primordial CMB to make use of the extra information. To examine our hypothesis, 
however, we can look at extreme cases. From equation ( l39l) . the largest deviations are 
found in models where e is large, and e and r] have opposite signs. For example, consider 
e = 0.015 and rj = —0.02, corresponding to spectral parameters r = 0.23, Us = 0.9 and 
a = 0.0. Figure [3] shows the Figures of Merit for the spectral and slow roll formulations 
for this parameter set, and the two values have begun to diverge once imax > 1500. 
This is a gedanken experiment, given both that we assumed foreground subtraction and 
that Ug — l and r are unrealistically large. This illustrates that second order consistency 
conditions are probably not testable by CMB data alone. However, it also gives cause 
for optimism that this lever arm effect can be exploited by combinations of CMB surveys 
and datasets sensitive to the primordial spectrum at shorter scales. 

Adding further decades in information from the power spectrum will increase 
the difference between the figures or merit. However, we also see that the running 
induced solely by e and t] being non-zero is proportional to e. Thus in Low-e models 
(corresponding to those with a negligible value of r), increasing the lever arm in In/c/fco 
will not give any extra advantage to slow roll reconstruction, relative to the spectral 
variables, unless we are lucky enough to see scale dependence ^ induced by a non-zero 
rj. Moreover, if we do observe a running in conjunction with a low value of r, then we 
will obtain an unambiguous measurement of ^. 
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Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal % 


Plk Abs 


Plk % 




0.024 


0.0000607 


0.253 


0.0000772 


0.322 


0.00018 


0.75 


^cam 


0.12 


0.000293 


0.244 


0.000503 


0.419 


0.0014 


1.17 


Ho 


72.0 


0.132 


0.184 


0.225 


0.312 


0.664 


0.922 


T 


0.164 


0.0026 


1.59 


0.00306 


1.87 


0.000637 


3.88 




0.9 


0.00437 


0.486 


0.00547 


0.607 


0.0116 


1.19 


e 


0.01 


0.000122 


1.22 


0.00021 


2.1 


0.00181 


18.1 


V 


0.0 


0.00151 




0.00206 




0.00331 




e 


0.0 


0.00174 




0.00216 




0.00285 




e 


0.01 


0.000122 


1.22 


0.000209 


2.09 


0.00179 


17.9 


V 


0.0 


0.000871 




0.00115 




0.00417 




e 


0.01 


0.000116 


1.16 


0.000193 


1.93 


0.000914 


9.14 



Table 3. Errors for model: e — 0.01, r/ = 0.0, ^ = 0.0. This model gives 49.5 e-folds 
of inflation after the fiducial scale leaves the horizon. Here C.V. refers to a cosmic 
variance limited survey. Ideal is the "straw man" satellite proposal of [52] and Plk 
refers to Planck (although this analysis assumes perfect foreground subtraction, and 
thus over-estimates Planck's capabilities, especially at small r). We give absolute and 
percentage errors for each parameter. In all cases we assume £max = 1500. The second 
block gives the forecast for a fit to {e,??,^}, the middle block for a fit to {e,??} and 
the bottom block give the results for a fit to e alone. As expected the constraints 
get tighter as the number of parameters is reduced. Finally the top block gives the 
forecasts for the other cosmological parameters, derived from the three parameter slow 
roU fit. 

6. Forecasts 

We use the derivative methodology outlined in [27] to compute the numerical derivatives 
using CAMB [26]. We consider the parameters 'class-by-class' as described in ^ 
beginning with the simplest possible case, the High-e 1-Parameter model. In each case 
where one or more of the parameters is zero, we consider both situations in which we 
fit for the parameter at its zero value and the alternative in which it is not included in 
the parameter set. This second case reduces the dimensionality of the parameter space 
by one and naturally leads to tighter restrictions on the other parameters. For all but 
the High-e 1-Parameter models, we work with a fiducial point chosen so that Ug = 0.97 
at the pivot. That is, r sets e while a sets ^ and then both e and ^ are used to set r] 
keeping Ug = 0.97. 

High-e 1-Parameter Models: This is perhaps the simplest possible model of inflation 
and effectively describes an m^0^ potential. However, the actual portion of the potential 
sampled by inflation is not specified, and if this turns out to be far from the minimum 
inflation must end via a hybrid transition. We consider e = 0.01 and e = 0.0025. Via 
equation IA.2I the first case has a further 49.5 e-folds of inflation after the fiducial scale 
has left the horizon and gives a tilt of = 0.96. The second case has 199.5 further 
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Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal % 


Plk Abs 


Plk % 


e 


0.0025 


0.0000522 


2.09 


0.000112 


4.48 


0.000826 


33. 


V 


0.0 


0.00151 




0.00206 




0.00369 






0.0 


0.00172 




0.00214 




0.00332 




e 


0.0025 


0.0000522 


2.09 


0.000112 


4.48 


0.000825 


33 


V 


0.0 


0.00083 




0.00108 




0.00266 




€ 


0.0025 


0.0000517 


2.07 


0.000109 


4.38 


0.000665 


26.6 



Table 4. Errors for model: e = 0.0025, rj = 0.0, ^ = 0.0. This model gives rise to 
199.5 e-folds of inflation after the fiducial scale leaves the horizon, and we have used 
the same conventions as Table O 



e-folds and a tilt of Ug = 0.99. This model has no obvious analogue in spectral variables, 
so we only present forecasts for the slow roll variables. As can be seen from Tables [3] 
and m this one parameter model can be very tightly constrained. However, the precise 
measurement of r is facilitated by the correlation of r and Ug, and the latter is very 
tightly constrained by measurements of the temperature anisotropics. 

Low-e 1-Parameter Models: In this case, the tensor signal is tiny, and |?7| ^ 0. Table E] 
shows the forecast errors for a model rj so that Us = 0.97 and a = 0, so ?7 = —0.015 
and ^ = 0.0. This scenario model leads to ~ 600 e-folds of inflation, so we would 
clearly require a hybrid transition. Note that when ^ and a are marginalized over, the 
relative constraints on — 1 and t] are identical. We make no forecast for e - this fit 
is analogous to that of a pure ACDM cosmology. This forecast includes the B-mode 
contribution from the lensed E-mode. We present the same forecast in Table [6], with 
the S-mode spectrum omitted from our analysis - that is all derivatives of C^^ are 
assumed to vanish. We see that in the absence of primordial tensors, a highly accurate 
5-mode measurement would significantly improve constraints on ujcdm, since the lensing 
is induced by the dark matter potential wells - although whether this level of accuracy 
is achievable in practice is of course a very different matter. 

We give forecasts for both slow roll and spectral variables. We construct the 
analogous spectral variables by computing them at the pivot. The resulting spectra 
are subtly different, since a two parameter slow roll model has a non-zero running if e 
is not minuscule; this effect will necessarily be absent in the spectral case. If we assume 
that the underlying cosmology is a Low-e 1-Parameter model, we are still free to include 
^ and a in our parameter set. In this case, a ~ —2^ and we see that the constraints on 
^ indeed match those on a. However, adding weakens the constrain on rj since Ug has 
a partial degeneracy between a and rj. This amounts to a rotation in our parameter 
space, and does not imply that slow roll constraints are intrinsically less accurate than 
the spectral constraints. We implement "Low-e" scenarios by using the full slow-roll 
hierarchy with e set to a very low value (10^^°). The results are insensitive to the actual 
value of e, provided it is far below the threshold of detectability. 
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Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal % 


Plk Abs 


Plk % 




0.024 


0.0000614 


0.256 


0.0000778 


0.324 


0.000182 


0.758 


^cdm 


0.12 


0.000311 


0.259 


0.00051 


0.425 


0.00142 


1.19 


Ho 


72. 


0.141 


0.195 


0.229 


0.317 


0.676 


0.939 


T 


0.164 


0.00281 


1.71 


0.00328 


2. 


0.00647 


3.94 


A{ko) 


0.9 


0.00468 


0.52 


0.00577 


0.641 


0.0117 


1.3 


Us - I 


-0.03 


0.00166 


5.53 


0.00215 


7.16 


0.00439 


14.3 


a 


0. 


0.00344 




0.00428 




0.0066 




Us - I 


-0.03 


0.00166 


5.53 


0.00212 


7.07 


0.00427 


14.2 


V 


-0.015 


0.00308 


20.6 


0.00404 


26.9 


0.00632 


42.1 




0. 


0.00346 




0.00432 




0.00669 




V 


-0.015 


0.000828 


5.52 


0.00106 


7.05 


0.00213 


14.2 



Table 5. Errors for model e = 10-^°, 77 = -0.015 and ^ = 0.0, where the CMB 
spectra include the -B-mode contribute from the lensed E mode. The conventions here 
match those of Table |3l although we also give forecasts for the corresponding spectral 
parameters. 



Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal 


% Plk Abs 


Plk % 


UJb 


0.024 


0.000066 


0.275 


0.0000831 


0.346 


0.000192 


0.801 


^cdm 


0.12 


0.000634 


0.528 


0.000755 


0.63 


0.00155 


1.29 


Ho 


72. 


0.271 


0.376 


0.33 


0.459 


0.741 


1.03 


T 


0.164 


0.00297 


1.81 


0.00338 


2.06 


0.00642 


3.92 


A{h) 


0.9 


0.00574 


0.638 


0.00657 


0.73 


0.00117 


1.3 


Us - I 


-0.03 


0.00204 


6.8 


0.00248 


8.23 


0.00459 


15.3 


a 


0. 


0.00377 




0.00443 




0.00653 




Us-l 


-0.03 


0.00198 


6.6 


0.00237 


7.9 


0.00457 


15.2 


T) 


-0.015 


0.00363 


24.2 


0.00436 


29.1 


0.00632 


42.1 




0. 


0.0038 




0.00448 




0.00662 




V 


-0.015 


0.00099 


6.6 


0.00119 


7.9 


0.00228 


15.2 



Table 6. Errors for model e — 10^^", r/ — —0.015 and ^ = 0.0. In this case we did not 
include any B-mode information in our forecasts. The forecasts for our measurement of 
uJcdm are considerably looser than those found in Table [5] even though the primordial B 
mode is entirely absent, since the amplitude of the B mode is now entirely dependent on 
the depth of the dark matter potential wells, and we are not exploiting this information 
here. 
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Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal % 


Plk Abs 


Plk % 


^0 


0.024 


0000607 


0.254 


0000772 

KJ • KJ KJ KJ KJ 1 1 w 


0.322 


00018 

KJ • KJ KJ KJ -L Ky 


0.75 


^cam 


0.12 


0.000305 


0.254 


0.000516 


0.43 


0.00141 


1.18 


Ho 


72.0 


0.136 


0.189 


0.23 


0.319 


0.67 


0.93 


T 


0.164 


0.00271 


1.66 


0.00318 


1.94 


0.0064 


3.9 




0.9 


0.00451 


0.501 


0.00558 


0.62 


0.0116 


1.29 


r 


0.24 


0.0026 


1.08 


0.00425 


1.77 


0.351 


14.6 




-0.03 


0.00167 


5.57 


0.00217 


7.23 


0.00427 


14.2 


a 


0.0 


0.00355 


- 


0.0044 


- 


0.00675 


- 


r 


0.24 


0.0026 


1.08 


0.00425 


1.77 


0.0348 


14.5 


Us-l 


-0.03 


0.00167 


5.57 


0.00213 


7.1 


0.00425 


14.2 


t 


0.015 


0.000158 


1.06 


0.000259 


1.73 


0.00213 


14.3 


V 


0.015 


0.00154 


10.3 


0.00209 


13.9 


0.00425 


34.1 




0.0 


0.00176 




0.00218 




0.00337 




e 


0.015 


0.000157 


1.05 


0.000257 


1.72 


0.00269 


13.9 


V 


0.015 


0.000911 


6.08 


0.00122 


8.11 


0.00483 


32.2 



Table 7. Errors for model: e = 0.015, 77 = 0.015, ^ = 0.0 



High-e 2-Parameter Models: If we consider only scalar perturbations High-e 2- 
Parameter models have a degeneracy in the e — 77 plane, which is broken by the limits on 
e derived from B-mode observations, and/or the running of the spectral index. As we 
showed in Figure [H the constraints on the slow roll parameters will tighten more rapidly 
than the analogous ones on the spectral parameters when e is relatively large thanks 
to the lever arm effect. In the forecast presented in Table [7] we have set e = 0.015 
and chosen 77 so that Ug = 0.97. However, since our forecasts have been done with 
^mnr = 1500 we do uot See the benefit of the lever arm effect here. 



Low-e 2-Pammeter Models: We now consider the scenario in which the inflaton 
potential is specified by r] and ^. That is, r = 0, Ug ^ I and a 7^ 0. In this case, the 
slow roll parametrization is (almost) identical to the spectral parametrization. There 
will be terms proportional to ^\n{k/ko)^, but these are never large over the range of 
scales probed by the CMB. In computing the errors for these models, we do not include 
information from lensing of the E-mode polarization into B-mode polarization, so we 
have somewhat over-estimated the likely errors for Ucdm- However, both the spectral 
and slow roll variables will be equally affected. 

The parameter ^ dominates the dynamics very quickly when ^ > [69], and thus 
a < 0. In Appendix A we present an analytic treatment of the slow roll dynamics for 
the Low-e 2-Parameter model, from which we can compute the remaining number of 
e-folds. Fixing Ug = 0.97 and then decreasing a from zero leads to a dramatic decrease 
in Nf., as shown in Figure RTl For rj = —0.015, ^ = 0.0 inflation ends after ~ 600 e-folds. 
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Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal % 


Plk Abs 


Plk % 


Us - I 


-0.03 


0.002 


6.67 


0.00247 


8.23 


0.00459 


15.3 


a 


-0.004 


0.00378 


94.4 


0.0044 


111. 


0.00653 


163. 




-0.0165 


0.0019 


11.6 


0.00228 


13.8 


0.00343 


20.9 




0.002 


0.0019 


97.1 


0.00225 


112. 


0.00331 


165. 




Table 8. Errors for model: e = 10 


-1°, r/ ^ -0.0165, C = 0.002 






Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal % 


Plk Abs 


Plk % 


Us - I 


-0.03 


0.002 


6.67 


0.00236 


7.87 


0.00459 


15.3 


a 


-0.01 


0.00378 


37.8 


0.00444 


44.4 


0.00655 


65.5 


Tj 


-0.019 


0.0019 


10.2 


0.00228 


12.2 


0.00344 


18.5 


e 


0.005 


0.00195 


39 


0.00225 


45 


0.00332 


66.4 




Table 9. Errors for model: e = 10 


-10, rj = -0.019, C = 0.005 






Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal % 


Plk Abs 


Plk % 


ris-l 


-0.03 


0.002 


6.67 


0.00244 


8.13 


0.00457 


15.2 


a 


-0.06 


0.00381 


6.36 


.00449 


7.48 


0.00672 


11.2 


V 


-0.037 


0.00182 


4.91 


0.00222 


6.02 


0.00335 


9.09 




0.03 


0.00189 


6.31 


0.00223 


7.44 


0.00333 


11.1 



Table 10. Errors for model: e = lO^^", r/ = -0.037, ^ = 0.03 



while ^ = 10^^ gives ~ 260 e-folds and ^ = 0.001 yields only ~ 55 efolds, even though 
the corresponding value of a at the pivot is very small. At this point we are beginning 
to bite into the e-foldings window of equation ( l25l) . Almost any prior on the number of 
e-folds will lead to tight constraints on ^ [8J. This constraint has has no analog when the 
power spectrum is specified in terms of Ug, and a, and allows slow roll reconstruction 
to impose much tighter constraints than fits to Ug, and a. 

For this set of parameters, we find that the smallest detectable value of ^ (assuming 
a cosmic variance limited measurement of the CMB) is ^ ~ 0.002. Consequently, if ^ 
can be detected in via fits to CMB data, the end of inflation was 'in sight' as the 
fiducial modes left the horizon, and we do not need to posit a hybrid transition. Given 
the difficulties of foreground subtraction, a constraint this tight would almost certainly 
require a combined fit to several orthogonal datasets. Furthermore, as pointed out by 
[To], the extra parameters we have ignored (neutrino mass, dark matter equation of 
state), can have a non- negligible effect on estimates of primordial spectrum parameters 
and thus forecasts. We demonstrate our ability to constrain ^ in Tables ISlfTOl In each 
case we forecast the errors with t] and ^ chosen so that Ug = 0.97, and with increasing 
values of ^. 
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Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal % 


Plk Abs 


Plk % 


r 


0.24 


0.0026 


1.08 


0.00424 


1.77 


0.0351 


14.6 


Us - I 


-0.03 


0.00167 


5.57 


0.00216 


7.2 


0.00427 


14.2 


a 


-0.002 


0.00349 


175 


.00436 


218 


0.00669 


335 


6 


0.015 


0.000158 


1.05 


0.000258 


1.72 




1 1 

±0 . -L 


V 


0.14 


0.00154 


11.0 


0.00209 


14.9 


n nnzLSQ 


lA Q 


c 


0.001 


0.00176 


176 


0.00218 


218 


u.uuoou 


oou 




Table 11. 


Errors for model: e ^ 0.015, 77 = 0.014, i = 0.001. 


This model 


gives ^ 3 




e-foldings of inflation after the fiducial scale leaves the horizon 






Parameter 


Fiducial 


C.V. Abs 


C.V. % 


Ideal Abs. 


Ideal % 


Plk Abs 


Plk % 


r 


0.24 


0.00257 


1.08 


0.00422 


1.78 


0.0348 


14.6 




-0.03 


0.00167 


5.57 


0.00216 


7.2 


0.00427 


14.2 


a 


-0.02 


0.00357 


17.9 


.00442 


22.1 


0.00683 


34.2 


e 


0.015 


0.000158 


1.06 


0.00026 


1.73 


0.00213 


14.2 


ri 


0.0075 


0.00152 


20.4 


0.00207 


27.8 


0.00509 


68.3 


i 


0.01 


0.00176 


17.6 


0.00218 


21.8 


0.00337 


33.7 



Table 12. Errors for model: e = 0.015, -q = 0.0075, ^ = 0.01. This model gives only 
~ 12 e-foldings of inflation after the flducial scale leaves the horizon. 



High-e 3-Parameter Models: Finally, we consider the case in which the first three slow 
roll parameters are non-zero. In this case e can be large, and the slow roll variables can 
again evolve within the observable window. We present results for two representative 
cases in Tables [11] and [121 However, the principal constraining power of slow roll 
reconstruction, relative to the usual spectral formulation, is again its ability to include 
restrictions based on the duration of inflation. For example, the model in Table [T2] would 
manage a bare 12 e-folds of inflation. If we found these central values for {e, r^,^}, we 
would know this description of the inflaton potential was either incomplete or that the 
primordial spectrum was not generated by slow roll inflation. 

7. Conclusions 

We have considered the constraints one can place on the inflationary parameter 
space via slow roll reconstruction. In the course of this work we also analyzed the 
precision of different approximation schemes that can be used to compute the underlying 
spectrum. As noted in f7], the primordial spectrum is easily computed by numerically 
solving the mode equations fll3til4p . In this case slow roll reconstruction contains no 
approximations, as the slow roll hierarchy itself is solved analytically [16]. For almost 
all models of interest, approximate solutions will be more than adequate given the 
precision of current data. The principal concern is that different schemes for computing 
the spectrum will impose implicit priors on the duration of inflation, and these will effect 



Constraining Inflation 



26 



the computed likelihood contours if they are not properly understood. This issue was 
investigated in Refs. [T71 [18] , which implement slow roll reconstruction with an exact 
solver for the perturbation spectrum; they find that the sole difference in likelihood 
contours using current data can be attributed to implicit priors on the duration of 
inflation. However, since the spectrum can be computed numerically at little cost, 
using the exact solver has the advantage of removing a possible source of uncertainty. 

In addition, we introduce a scheme for classifying inflationary models, based on the 
order at which the slow roll hierarchy is truncated, and whether the e parameter is small 
enough for its contribution to the inflationary dynamics to be ignored. We refer to these 
models as "Low-e, N-parameter" and "High-e, N-parameter" models. In the former case, 
we can impose the slow roll prior while ignoring a possible tensor contribution to the 
CMB, in analogy with the usual ACDM parameter set of the concordance cosmology. 
Physically, e responds to the ratio of the amplitude of tensor and scalar perturbations, 
which reflects the energy scale at which inflation occurs. It is possible (and some might 
say likely, if one expects inflation to be consistent with the microphysical constraints of 
string theory) the inflationary scale is well below the GUT scale. In this case e and the 
amplitude of the tensor spectrum must be vanishingly small. In this limit e effectively 
decouples from the slow roll hierarchy. 

Slow roll reconstruction can lead to tighter constraints on the parameter space 
than those found with the usual spectral indices and runnings. In particular, slow roll 
reconstruction allows us to impose constraints on the total duration of inflation, because 
slow roll reconstruction requires one to compute the duration of the inflationary era for a 
given set of parameters. If this value is very low - as it can be for apparently reasonable 
parameters - we do not have a self-consistent description of inflation. This effect is 
clearly seen in implementations of slow roll reconstruction [6l [7]. Constraints based 
on the duration of inflation imply a sharp cut on the relevant parameter space, and 
this effect is not captured by our Fisher matrix forecasts which assume a Gaussian 
likelihood. These constraints simply encode our assumption that we are using enough 
slow-roll parameters to provide a full description of the inflationary potential, and do 
not introduce any significant new or ad hoc assumptions into the parameter estimates. 
The impact of this constraint could be forecast by running Markov Chains on simulated 
data for a proposed experiment. We sketch an improved "e-dependent" algorithm for 
imposing the "e-fold prior" that will be explored in detail in future work. 

Slow roll reconstruction can also lead to improved constraints via the lever-arm 
effect, which arises when the slow roll variables evolve appreciably in the range of scales 
open to cosmological experimentation. Higher order slow roll consistency relations 
are automatically incorporated into slow roll reconstruction, and the significance of 
these terms grows rapidly as once increases the range of wavelengths over which we 
have information about the primordial spectrum. Our error forecasts for the slow 
roll parameters show that it is very unlikely this effect can be exploited by datasets 
containing only CMB information. However, high redshift 21 cm data may probe the 
the primordial spectrum at scales unreachable via the CMB. We thus conjecture that 
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the combination of CMB and 21 cm data may test these consistency relations directly, 
fully realizing the power of slow roll reconstruction. This conjecture is supported by 
the recent analysis of [71], who project that 21 cm data could provide a 4(T detection of 
the intrinsic running predicted by simple inflationary models. This is the same running 
one expects from higher-order consistency conditions, which is the basis of the lever 
arm effect we described earlier. There is thus good cause for optimism that slow roll 
reconstruction will be able to exploit this information to significantly tighten constraints 
on the inflationary parameter space and we plan to explore this issue in detail in a 
future publication. If inflation occurs well below the GUT scale, e is very small, and 
the tensor modes are likely to be unobservable. The running in the spectral index 
induced by the scale dependence of the slow roll parameters is also proportional to e. 
Consequently, this effect is strongly correlated with the inflationary scale. Finally, BBO 
style direct detection experiments for gravitational waves probe the primordial spectrum 
at wavelengths some 30 e-folds smaller than those in the CMB [72| - and in this case 
the running of the spectral indices will be unmistakable if the tensor amplitude is large 
enough to be seen, and slow roll reconstruction will be well-placed to take advantage of 
these effects. 

If the tensor to scalar ratio is vanishingly small, the tilt of the scalar spectrum is 
fixed by rj, if |^| <^ jr]]. If r] = —0.015, corresponding to = 0.97, inflation lasts for 
much longer than 60 e-foldings. If our chosen set of variables provides a complete 
description of the spectrum, inflation must end via an some abrupt change in the 
potential, which cannot be encoded in the truncated slow roll expansion. Conversely, 
in models where e is non-trivial, or the ^ parameter is even slightly greater than zero, 
inflation will typically end within 60 e-folds, and often less. In this case, the inflaton 
can 'see' the end of inflation well before it occurs. Further, if e is constrained to be 
very small by a tight upper bound on r, any running detected in the spectrum must be 
induced by an ^ parameter, which is equivalent to learning that V"'{(f)) is non-trivial. 

In summary, this paper provides explores the theoretical basis of slow roll 
reconstruction, as originated by Easther and Peiris, and extends our understanding 
of the properties of the slow roll hierarchy. Slow roll reconstruction requires no explicit 
approximations in the treatment of the inflationary dynamics and exploits all possible 
sources of astrophysical information. Consequently it provides an optimal solution to the 
problem of reconstructing the inflationary potential from astrophysical data. Whether 
the universe will be kind enough to permit this program to be implemented in practice 
- and whether the generation of the primordial spectrum can in fact be attributed to 
an inflationary model which can be described in terms of a single, slowly-rolling fleld - 
remains to be seen. 
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Appendix A. Solutions to the Slow Roll Hierarchy 

High-e 1 -Parameter: The case where we truncate the flow equations at first order 
corresponds to a potential that is quadratic, at least around the pivot. We can solve 
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the flow equations exactly in this case, irrespective of the value of e: 

where eo is the value of e at the fiducial scale (where we have set = 0). We find that 
inflation ends after 

1-eo 



end 



2er 



(A.2) 



The value of e can be tuned to obtain the right amount of inflation to satisfy the 
matching condition, equation ( l25i) . However, as we have discussed above, in a one 
parameter model, e is tightly constrained by the tilt, tensor-scalar ratio and tensor tilt. 
Any tension between these quantities will result in the need for non-zero higher order 
HSR parameters. 

Low-e 1-Parameter: See section |3] for details of this case. 



Low-e 2-Parameter: In this situation, the third flow equation is 

In the limit e ^ 1 we can solve equations (l26l) and (1A.3I) to obtain 

1 



-^^{2B,-N) 



iv^(2i?. 



N) 



r^[N) = - A/2A^tan 
^(N) = A^tan^ 

where = — ^rj^ and = arctan(— r/o/ a/2v4^)/ >/2A^. Finally, 
e(A^) = C, sec 



(A.3) 

(A.4) 
(A.5) 



-V^(2i?,-Ar) 



where C, = e^/ sec^{^/2A^B^). We can now solve for where inflation ends (where e 
in this approximation: 

2 arccos icl^'^ 



1) 



- A^ 



y2A^ 



2 Br. 



(A.6) 



The upper panels of figure lAll show the quality of the approximation, which only gets 
better as one decreases the value of e. The lower panels of figure lAll shows the number 
of e-folds as a function of the slow roll parameters rj and ^ and the tilt and running of 
the scalar spectral index respectively. The end of inflation in these models is completely 
specified by the initial values of the slow roll parameters. 

Using a value of ^ that gives a running anywhere near the centroid of the WMAP 
data [51j leads to a short period of inflation. If the current weak evidence for running 
detected in the WMAP data persists, higher order slow roll parameters are required to 
produce enough inflation. This is consistent with the findings of [8]. 
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Figure Al. Slow roll hierarchy truncated at third order (^Xh = 0) and ^ = 0.001. 
The left panel shows the evolution of e while the right panel show rj; the red curve is 
the exact result, the blue is the approximation. The lower left plot shows the amount 
of inflation as a function of the parameters ry and ^, while the lower right plot shows 
the same quantity as a function of the spectral indices. 




Figure A2. Slow roll hierarchy truncated at fourth order ^Ah = 0. At the pivot, 
T] = -0.02, f = 0.03, ^Xh ^ -0.0001 and e = lO"'^. Left panel show the evolution of 
e, the red curve is the exact result, the blue is the approximation. 



Low-e 3-Parameter If we include ^Xh there is no obvious analytic solution. However, 
if ^Xh is small, and noting that rj is also small and remains small throughout inflation, 
we employ a further approximation: 

In this limit rj obeys Airy's equation, with solution: 



viN) 



h\ / CiAi-[x(N)]+Bi-[x(N)] \ 
J { CiAi [x(N)] + Bi [x(N)] ; ' ^ ■ > 
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where 



x{N) = T 



±3A 



H 



1/3 



A, 



(A.9) 



and Ci is a constant given by: 

/ 

Ci = - 



Bi' 



+ 



-1/3 



VoBi [x(0)] 



Ai' 



X 



+ 



-1/3 



r/oAi [x 



(A.IO) 



Now we can also solve for e(A^) to obtain 

e{N) = — J (A.ll) 

(CiAi[x(N)]+Bi[x(N)]f 

where C2 is given by: 

C2 = eo (CiAi [x(0)] + Bi [x(0)])' (A.12) 

Figure IA2I shows the quality of the approximation, even for large values of the running, 
or ^. The value of ^ here is taken to coincide with a value of a = —0.06, i.e. ^ = 0.03. 
Notice how quickly inflation ends; barely 10 e-folds of inflation are achieved. 

These expressions are of limited use in practice. If ^Xh is large enough to prevent 
^ from dominating the dynamics, it must be of the same magnitude as ^, and equation 
( ]A.7|) is no longer a good approximation. The effect of including ^Xh depend on its sign. 
If '^Xh < 0, inflation ends more rapidly, for fixed {e, ^}. If < ^Xh ^ 1, we see a slight 
increase in the number of e-folds. As one increases ^Xh, one encounters a critical point 
is reached where e has a maximal value less than unity and then begins to decrease. 



